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RESUMEN 

Aunque la materia barotropica no constituye un modelo realista para estrellas magneticas, seri'a interesante 
confirmar una conjetura reciente que establece que las estrellas magneticas con ecuacion de estado barotropica, 



serian dinamicamente inestables (Reisenegger 2009). En este trabajo construimos un conjunto de equilibrios 
barotropicos, los cuales pueden ser finalmente testeados usando un criterio de estabilidad. Una descripcion 
general de las ecuaciones de MHD ideal que gobiernan estos equilibrios es revisada, permitiendo tanto una 
componente poloidal, come una componente toroidal del campo magnetico. Un nuevo codigo numerico en 
diferencia finita es desarrollado para resolver la llamada ecuacion de Grad-Shafranov que describe el equilibrio 
de estas configuraciones, y algunas propiedades de los equilibrios obtenidos son brevemente discutidas. 

ABSTRACT 

Although barotropic matter does not constitute a realistic model for magnetic stars, it would be interesting 
to confirm a recent conjecture that states that magnetized stars with a barotropic equation of state would be 



dynamically unstable ( Reisenegger 2009 ) . In this work we construct a set of barotropic equilibria, which can 
eventually be tested using a stability criterion. A general description of the ideal MHD equations governing 
these equilibria is summarized, allowing for both poloidal and toroidal magnetic field components. A new 
finite-difference numerical code is developed in order to solve the so-called Grad-Shafranov equation describing 
the equilibrium of these configurations, and some properties of the equilibria obtained are briefly discussed. 
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1. BARO... WHAT? 

Barotropic equations of state, where pressure 
is a function solely of density, are often assumed 
to describe the matter within magnetic stars in 
ideal magnetohydrodynamic (MHD) equilibrium 



(Yoshida & Eriguchi 


Lander & Jones 


2009) 



2. MHD EQUILIBRIA: THE 
GRAD-SHAFRANOV EQUATION 

In the ideal MHD approximation, a magnetic star 
may be considered as a perfectly conducting fluid in 
dynamical equilibrium described by the Euler equa- 
tion. 



the range of possible equilibrium configurations, 
and does not strictly represent the realistic stably 
stratified matter within these objects, which is 
likely to be an essential ingredient in the stability 
of magnetic fields in stars (Reisenegger [2009 ). 
With this in mind, it is interesting to carry out 
the pedagogical exercise of checking whether the 
unrealistic barotropic equilibria are really stable or 
not. This work is focused on obtaining a wide range 
of these equilibria and study their main properties, 
as a starting point to study their stability. 
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V_P + pV$ = - J X B, 



(1) 



where the right-side is the Lorentz force per unit 
volume. Considered objects have a very large fluid 
pressure P {P ^ GM^ /R^, M being the mass and 
R the radius), to magnetic pressure B'^/Stt ratio 
{B being a char acteristic magneti c field strength), 
SttP/S^ - 10** (|Reiseneggerl|2009|), which suggests 



that magnetic forces may be balanced by a slight per- 
turbation of an unmagnetized spherical background 
equilibrium. This implies that, as an approximation, 
we can consider the star as spherical with negligible 
deformations due to the magnetic forces. In addi- 
tion, if axial symmetry is assumed, and spherical co- 
ordinates (r, 9, (j)) are used to describe the model, all 
scalar quantities are independent of the azimuthal 
coordinate, and the magnetic field may be expressed 
as the sum of a poloidal (meridional field lines) com- 
ponent, and a toroidal (azimuthal field lines) com- 
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Fig. 1. Meridional cut of a star bearing an axisymmetric 
magnetic field. The bold curve is the surface of the star, 
while the thinner curves are poloidal field lines. The 
toroidal component of the magnetic field may lie only in 
regions where the poloidal field lines close inside the star 
(gray region). 



ponent, each determined by a single scalar function, 



B = B 



pol ' 



Bt 



Va(r,6') X V0 + /3(r,6l)V</>, (2) 



which turn out to be constant along their respective 



field lines ( Chandrasekhar & Prendergast 1956). Un- 



der this symmetry, the azimuthal component of the 
magnetic force per unit volume must vanish, which 
implies a functional relation between these scalar 
functions, f3{r,9) = /3{a{r,9)). In this way, both 
a and /? are constant along field lines and, if vacuum 
is assumed outside the star, the toroidal field may 
lie only in regions where the poloidal field lines close 
within the star (Figure [ij . On the other hand, if a 
barotropic equation of state, P = P{p), is assumed, 
the Lorentz force per unit mass must be the gradient 
of some arbitrary function x(r, 9), which turns out to 
be a function of a as well, x(?', Q) — X {(^{^^ Using 
all this formalism, a non-linear elliptic partial differ- 
ential equation is found to be the master equation 
governing the equilibrium of a barotropic MHD equi- 
librium, the so-called Grad-Shafranov (GS) equa- 
tion. 



Q^2 



1 da 
"2 d9 \sin9d9 



■ PP' 4- sin^ 6'px' = 
(3) 



( |Grad fc Rubin|1958[ |Shafranov|1966l ) where ] crimes 
stand for derivative with respect to the argument, 
and both j3 = 13(a) and x — x(q^) ^re two arbi- 
trary functions, whose form may be chosen depend- 
ing on the particular magnetic configuration of in- 
terest. Under the assumption of weak magnetic field 
discussed in ^JT] the density p appearing in the GS 



equation may be replaced by its non-magnetic back- 
ground counterpart, p = p{r), such that we solve for 
the magnetic functions for a given density profile, 
instead of considering the more difficult problem of 
solving self-consistently for the magnetic functions 
and for the fluid quantities. 

3. NUMERICAL SOLUTIONS 

Outside the star, a corresponds to an infinite su- 
perposition of multipoles, which is the general so- 
lution of the GS equation with both (3 — and 
p — 0. We have implemented a finite-difference code 
to solve numerically the GS equation inside the star, 
for arbitrary choices of f3{a) and x(q:)- Solutions 
are matched to the exterior expansion by demand- 
ing continuity of a and its derivatives (related to 
the magnetic field components), in order to avoid 
surface currents. After testing our code, we ob- 
tained barotropic equilibria for the particular case 
x(a) = a, p{r) = pc{l — r'^ / R^) and 



/?(«) 



s(a - asY 




as <a 
a < a,, 



(4) 



where s is a free parameter accounting for the rela- 
tive strength between the poloidal and the toroidal 
component. In the definition below, an exponent 
larger than 1 was chosen in order to prevent a 
discontinuous /?' at the layer between regions with 
and without toroidal field; it is found that larger 
exponents than 1.1 give a smaller toroidal field 
strength. Also, = «(_/?, 7r/2) stands for the value 
of a along the largest poloidal field line closing 
within the star, being R the stellar radius, so the 
toroidal field lies in the region inside the curve 
as = a{r, 9) only. Figure [2] shows some numerical 
results for the particular j3{a) in Eq. Q. Black 
lines correspond to poloidal field lines with 0.2q!s, 
0.4as, 0.6as, 0.8as, l.Oag, l.OSag and l.lSas, 
respectively, whereas the color map accounts for the 
strength of the toroidal field. In turn. Figures |3]|4] 
show the strength of the magnetic field along the 
axis and the equatorial line for these equilibria. 



4. DISCUSSION 

All the equilibria found consist of a mixed 
poloidal-toroidal field with a dominant poloidal 
component in the magnetic energy -Emag- For the 
cases studied so far, the energy stored in the toroidal 
component -Etor is only a few percent of the total 
magnetic energy, even in cases where the maximum 
strength of the toroidal field is comparable to that 
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Fig. 2. Numerical equilibria found with our code. Left: 
s = 10, with BtorZ-Bmag ~ 0.5%. Right: s = 35, with 

-Etor/£^mag ~ 3.2%. 
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Fig. 3. Magnitude of the magnetic field for the equilibria 
shown in Figure [2] with s = 10. The maximum toroidal 
strength is about one order of magnitude smaller than 
the poloidal one. 

of the poloidal component: the larger the toroidal 
field, the smaller the volume where it lies. This 
small contribution to the energy has already been 
reported in the literature, but assuming a purely 
dipolar magnetic field outside the star (Lander 
& Jones 2012). Our code, allowing an arbitrary 



number of multipoles, seems to indicate that higher 
multipoles do not contribute significantly to the 
energy of these equilibria, at least not for small 
to moderate values of s. It is desirable to study 
this fact in more detail and confirm, for instance, 
whether a global maximum for i^torZ-Emag exists, 
already reported using general-relativistic MHD 
( Ciolfi et al.||2009 |. Once we obtain a wide range of 
relevant equilibria with consistent physical choices 
of the arbitrary magnetic functions, their dynamical 
stability may be analyzed using either a perturbative 
analysis or numerically solving the time-evolution 
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Fig. 4. Magnitude of the magnetic field for the equilibria 
shown in Figure [2] with s — 35. Both the poloidal and 
the toroidal maximum strength are of the same order of 
magnitude. 

of such configurations. 
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